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TOO NAIVELY IMPLEMENTED CONFINEMENT MIGHT 
ALLOW FOR UNSTABLE BOUND STATES 

The perhaps most popular three-dimensional reduction of the Bethe-Salpeter framework 
for the description of bound states within quantum field theories is the Salpeter equation, 
found in the instantaneous limit of the Bethe-Salpeter formalism if assuming in addition 
free propagation of the bound-state constituents. Unfortunately, depending on the chosen 
Dirac structure of the Bethe-Salpeter kernel (which encodes all interactions between the 
bound-state constituents) even for Salpeter equations with confining interactions, arising, 
e.g., from quantum chromodynamics, the predicted — supposedly stable — bound states 
exhibit instabilities, probably of similar nature as those observed in Klein's paradox. The 
observation of this kind of instability has been reported by numerous but to a large extent 
numerical studies investigating confining interaction potentials of in configuration space 
either strictly linear shape or a form interpolating between harmonic-oscillator and linear 
behaviour [1,2]. Generally, for linearly confining Bethe-Salpeter kernels being mixtures 
of Lorentz scalar plus time-component Lorentz vector, stability can be assured only if the 
Dirac structure of this kernel is predominantly of time-component Lorentz- vector nature. 

Building on experience acquired in previous analyses, focused to the simpler "reduced 
Salpeter equation" [3, 4] and its improvement [5, 6] by including dressed propagators for 
all bound-state constituents [7-9], this investigation aims at rigorous analytic proofs [10] 
of the stability of solutions of the (now) full Salpeter equation with confining interactions 
(such as harmonic-oscillator potentials) of frequently used Lorentz structures, in order to 
identify all those kernels for which bound-state stability can be taken as granted ab initio. 
We regard bound states as stable if their energy (or, in their center-of-momentum system, 
mass) eigenvalues belong to a real, discrete portion of the corresponding spectrum that is 
bounded from below. Such discussion might provide further insight into the reasons why, 
for Lorentz structures different from a time-component Lorentz vector, instabilities arise. 



FULL (IN CONTRAST TO REDUCED) SALPETER FORMALISM 
FOR RELATIVISTIC FERMION-ANTIFERMION BOUND STATES 



Assuming, as usual, the Lorentz structures of the effective couplings of both fermion and 
antifermion to be represented by identical Dirac matrices (generically called T hereafter) 
and denoting the associated Lorentz- scalar interaction function by Vr(p, q) , our Salpeter 
eigenvalue equation governing a chosen fermion-antifermion bound state of mass M and 
distribution 4>(p) of internal momenta p (its "Salpeter amplitude") reads in the rest frame 

* (p) = / I g jEFr(p , q) (^(P)^*(q)rA-(p) ft 



M-2E(p) 
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with one-particle kinetic energy E (p) and energy projection operators A ± (p) , defined by 



E(p) = \/p 2 + m 2 , P=|p|, and A ± (p) 



E(p)±Yo(y-p + m) 



2E(p) 

here m labels the common mass of the bound fermion and the associated antiparticle. The 
projector structure of Eq. (1) constrains all its solutions, the Salpeter amplitudes 4>(p), to 

A ± (p)4>(p)A ± (-p)=0. (2) 



GENERAL FEATURES OF THE EIGENVALUE SPECTRA OF 
FULL SALPETER EQUATIONS 

The structural equivalence of Salpeter's equation (1) to the random-phase-approximation 
equation [11], familiar from the investigation of collective excitations in nuclear physics, 
or direct inspection allow one to identify various characteristics common to all solutions: 

• The random-phase-approximation structure of the Salpeter equation guarantees that 
the squares M 2 of the mass eigenvalues are real. In general, the spectrum itself is not 
necessarily real and, even where it is proven to be real, it is not bounded from below. 

• The most important applications of the instantaneous Bethe-Salpeter formalism are 
those which adopt interaction kernels only composed of momentum- space potential 
functions Vr(p, q) and Dirac couplings T that satisfy (p, q) = Vr(p, q) = Vr(q, p) 
and 70 fo = ±r. Their spectra of mass eigenvalues M in the complex-M plane are 
just unions of real opposite-sign pairs (M, —M) and/or imaginary points M = —M* . 

Since eigenvalues embedded in a continuous part of the spectrum may cause instabilities, 
the nature of the entire spectrum proves to be crucial: Apart from demanding eigenvalues 
to be real and bounded from below the bound-state stability sought is established if either 
eigenvalues and continuous spectrum are disjoint or the spectrum is purely discrete at all. 

As a consequence of the constraint (2), any Salpeter amplitude 4>(p) may be expanded 
in terms of at most eight independent components called, say, 0;(p) , i = 1 , . . . , 8. We start 
this quest for bound-state stability at a system requiring the least number of components. 



HARMONIC-OSCILLATOR CONFINEMENT SIMPLIFIES THE 
SALPETER INTEGRAL EQUATION TO AN EASIER-TO-HANDLE 
SYSTEM OF RADIAL EIGENVALUE DIFFERENTIAL EQUATIONS 



All the instabilities under consideration are expected to show up first in the pseudoscalar 
sector [1]. Consequently, the primary targets of all analyses of the present kind are bound 
states with spin-parity-charge conjugation assignment J PC = h . A Salpeter amplitude 
3>(p) describing such a state involves just two independent components (p) and 02 (p) : 



*(p) 
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is the unique form of any Salpeter amplitude 4>(p) for a fermion-antifermion bound state 
of total spin 7, parity P= ( — l) y+1 , and charge-conjugation quantum numberC = (— \) J . 

Let our interaction kernel be of convolution type [Vr(p,q) = Vr(p — q)L arising from a 
central potential V(r), where r = |x|, in configuration space. Using a harmonic-oscillator 
potential V(r) =ar 2 , where a ^ avoids triviality, to exemplify our line of reasoning, the 
Salpeter equation (1) reduces to a set of second-order differential equations utilizing only 
a single differential operator (which is just the Laplacian A = V • V, acting on £ = states) 



D = 



d 2 2d 
dp 2 p dp 



Assuming massless bound-state constituents, m = 0, facilitates our study of this problem. 

In our analysis interaction kernels of time-component Lorentz-vector, = y° <g> y°, 
and Lorentz-scalar, r<g>r = 1 <g) 1, nature may be discussed simultaneously if introducing 

+ 1 for T ® T = y° cg> y° (time-component Lorentz-vector interactions) 

— 1 for r eg) r = 1 <g) 1 (Lorentz-scalar interactions) 



a = 



as a discriminating parameter. By factorizing off all dependence on angular variables, the 
Salpeter equation for harmonic-oscillator interactions of time-component Lorentz-vector 
or Lorentz-scalar Dirac structure reduces to a system of two radial differential equations: 



(2p-aoD)^2(p)=M^(p) 



Mp) = M <h(p) 



2p-a[D-jp 



(3) 



Adding the sets for a = + 1 and o = — 1 , we find the Salpeter equation for Lorentz-scalar 
plus time-component Lorentz-vector mixing, r®r = ^y°®y° + 77 1 ® 1 , where t, , r\ 6l: 



[2p-a^-J])D]Up)=MUp) 
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p< 



MP) =M<h(p) . 



(4) 



The units of the momentum p can be chosen, without loss of generality, such that \a\ = 1. 

However, in spite of the considerable technical simplification achieved by reduction of 
the Salpeter integral equation to a set of differential equations, for the kernels of different 
Lorentz structure we shall perform the spectral analyses still due on a case-by-case basis. 



CONFINING INTERACTION KERNELS OF TIME-COMPONENT 
LORENTZ-VECTOR DIRAC STRUCTURE (T <g> T = f <g> f) 



For the time-component Lorentz- vector Dirac couplings (c = +1), the mere existence of 
bound states requires the coupling a to be positive: a > 0. In this case, an entirely analytic 
proof of the (so far basically numerically established) stability of the bound states may be 
constructed. Expressed by the positive self-adjoint operators on the Hilbert space L 2 (M 3 ) 

A = -A + 2r = A f >0, B = -A + 2r+ 4 = # + > , r = Ixl , 

our radial Salpeter equation (3) with harmonic-oscillator confining interactions is clearly 
equivalent to the (vanishing angular momentum sector of the) matrix eigenvalue problem 



Reiteration of the L.H.S. operation yields an equivalent problem for real eigenvalues M 2 : 

\B \( 2 ( fi\ A5/ 1= M 2 /i, 

BA ){ f 2 ) \h) BAf 2 =M 2 f 2 . 

The spectral theorem for self-adjoint operators allows one to define the unique positive 
self-adjoint square root of A: A 1 / 2 = (A 1 / 2 )^ > 0. Settings = A 1 / 2 f 2 converts the relation 
BA f 2 = M 2 f 2 to an eigenvalue equation, Qg = M 2 g, of the positive self-adjoint operator 
Q = A l l 2 BA l l 2 = > 0. Obviously, our operators A and B satisfy the inequality A < B; 
left and right multiplication by A 1 / 2 transforms this inequality into A 2 < A 1 / 2 5A 1 / 2 = Q. 
A theorem about the spectrum of a Hamiltonian with potential increasing beyond bounds 
claims that a Schrodinger operator H = —A + V, defined as sum of quadratic forms, with 
locally bounded, positive, infinitely rising potential V [V (x) — y +<*> for r — > °°] has purely 
discrete spectrum; application of this theorem guarantees that the spectrum of A is purely 
discrete. The spectral theorem for self-adjoint operators allows for a representation of A 2 
in terms of the same projection- valued spectral measure as A; consequently, the spectrum 
of A 2 is also entirely discrete. Trivially, the positive operators A 2 and Q are bounded from 
below. Combining [12-17] the characterization of all discrete eigenvalues of an arbitrary 
self-adjoint operator bounded from below by the minimum-maximum principle with the 
operator inequality A 2 < Q shows that the spectrum of our real squared mass eigenvalues 
M 2 of Q must be purely discrete and therefore the spectrum of bound-state masses M too. 

This stability of/w//-Salpeter bound states, proven for time-component Lorentz- vector 
harmonic-oscillator interaction with positive strength a > 0, excludes the possibility that, 
in the case a = +1, instabilities of bound states are induced by energy eigenvalues being 
embedded in a continuous spectrum of the Salpeter operator controlling the bound states, 
instead of belonging to its discrete spectrum, as expected for true confinement. The proof 
cannot be transferred to a < as in this case the counterparts of our operators A and B are 

2 

A = A + 2r, 5 = A + 2r-^. 
For both operators A and B, positivity is lost. This fact invalidates most steps of our proof. 



CONFINING INTERACTION KERNELS OF LORENTZ-SCALAR 

(r <g> r = 1 <g> 1) OR LORENTZ-PSEUDOSCALAR (T <g> T = y 5 <g> y 5 ) 

DIRAC STRUCTURE 



In the limit of the bound- state constituents being exactly massless, the Salpeter equations 
for Lorentz- scalar [i.e., a = — 1 inEq. (3)] and Lorentz-pseudoscalar interaction kernels 
(must) become identical, as the anticipated manifestation of chiral symmetry. The matrix 
eigenvalue problem is formulated in terms of the operator pairs (A, B) for a > and (A, B) 
for a < 0. In both cases, one of these operators is not positive, which spoils our reasoning. 



LINEAR COMBINATIONS OF BOTH TIME-COMPONENT 
LORENTZ-VECTOR (r<g> T = f <g> y°) AND LORENTZ-SCALAR 
(r<g>r = 1 <g> 1) CONFINING INTERACTION KERNELS 

Introducing the abbreviations a = a £ e K, /3 = a rj G K, the two operators entering in the 
Salpeter equation (4) for some mixture of scalar and time-component vector kernels read 

sf= (a-j8)(-A)+2r , 38= (a + j8) ^-A+^j +2r , r = |x| . 

Positivity of £o?/z (symmetric) operators ^ and presence of both derivatives demands 
a — /3 > and a + /3 > 0. These two relations restrain the couplings a and /3 to the range 

a = |a|>|j8|>0 -l<-<+l and a>0. (5) 

A proof similar to the above for a pure time-component Lorentz- vector kernel establishes 
stability irrespective of the relative sign of the two contributions of unequal Lorentz type: 
the dominance of the time-component Lorentz- vector kernel guarantees stability. Table 1 
illustrates these findings for various parameter ratios /3/a within the tolerable region (5). 

TABLE 1. Lowest-lying positive mass eigenvalues (in units of tya) of the Salpeter equation 
(4) for several values j3 / a of our couplings in the stability-compatible range — 1 < j3/« < +1. 
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6.15 
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7.07 


7.16 


7.15 
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6.87 


6.49 


4.68 


2 


7.63 


8.78 


9.14 


9.30 


9.33 


9.24 


9.01 


8.55 


6.14 


3 


9.01 


10.54 


11.02 


11.23 


11.28 


11.18 


10.92 


10.38 


7.45 
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10.31 


12.18 


12.75 


13.01 


13.07 


12.97 


12.68 


12.06 


8.65 


5 


11.55 


13.73 


14.38 


14.67 


14.75 


14.65 


14.32 


13.63 


9.77 


6 


12.74 


15.19 


15.92 


16.25 


16.35 


16.23 


15.87 


15.11 


10.83 


7 


13.88 


16.59 


17.39 


17.76 


17.86 


17.74 


17.35 


16.52 


11.84 
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14.98 


17.93 


18.81 


19.21 


19.33 


19.19 


18.77 


17.87 


12.81 


9 


16.04 


19.23 


20.19 


20.64 


20.78 


20.63 


20.16 


19.18 


13.74 



SUMMARY, FINDINGS, CONCLUSIONS, AND PERSPECTIVES 



Salpeter equations with potential functions rising to infinity in configuration space do not 
automatically predict stable bound states: for this to happen, the Lorentz behaviour of the 
involved Bethe-Salpeter kernels too is crucial. Truly confining interaction kernels can be 
singled out rather systematically by requiring the emerging bound-state energy spectra to 
be both real and discrete. This task is comparatively easy if the Salpeter equation reduces 
to a differential equation (as for the harmonic oscillator) but might be conducted for more 
general (say, linearly confining) potentials if all integral-equation issues can be mastered. 
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